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SUPERCONDUCTIVITY IN THE "HOT SPOTS" MODEL 
OF THE PSEUDOGAP STATE. 

E.Z.Kuchinskii, N.A.Kuleeva, M.V.Sadovskii 

Institute for Electrophysics, Russian Academy of Sciences, Ural Branch, Ekaterinburg 620016, Russia. 

Abstract — We analyze the anomalies of superconducting state in the model of pseudogap state induced by fluctu- 
ations of short - range order of "dielectric" (AFM (SDW) or CDW) type, and based on the scenario of "hot spots" 
formation on the Fermi surface, with the account of all Feynman graphs for electron interaction with pseudogap 
fluctuations, leading to strong scattering around the "hot spots". We determine the dependence of superconducting 
critical temperature T c on the effective width of the pseudogap, correlation length of short - range order and con- 
centration of nonmagnetic impurities. We also discuss possible connection of these results with the general form of 
the phase diagram of superconducting cuprates. 



One of the main problems in the physics of high - 
temperature superconducting cuprates remains the 
theoretical explanation of the typical phase diagram 
of these compounds [1]. Especially important is the 
clarification of the nature of the pseudogap state, 
observed in a wide region of temperatures and car- 
rier concentrations [2], which is obviously crucial 
for the understanding of electronic properties both 
in normal and superconducting states. Despite the 
continuing discussions, the preferable "scenario" of 
the pseudogap formation seems to be based on the 
model of the strong scattering of electrons by an- 
tiferromagnetic (AFM, SDW) short - range order 
spin fluctuations [2,3]. 

In a recent paper [6] we have presented micro- 
scopic derivation of Ginzburg - Landau expansion 
and studied the influence of these pseudogap fluctu- 
ation on basic superconducting properties (for both 
s— and d-wave pairing) in the model of "hot spots" 
the Fermi surface. Similar analysis using Gorkov's 
equations was given in Ref. [7]. 

In the model of nearly antiferromagnetic Fermi 
- liquid [4] electron interaction with spin fluctua- 
tions is usually described by dynamic susceptibility, 
characterized by correlation length £ and frequency 
u> s f of spin fluctuations, which are to be determined 
from experiment and can depend both on carrier 
concentration and temperature. Both dynamic sus- 
ceptibility and effective interaction are maximal (in 
momentum space representation) in the vicinity of 
vector Q = (7r/a, n/a) (a - lattice constant), which 
leads to the appearance of "two types" of quasipar- 
ticles — "hot" one with the momenta close to the 
points on the Fermi surface, connected by scattering 
vector ~ Q, and "cold" one with the momenta close 
to the parts of this surface surrounding diagonals of 
the Brillouin zone [2,4,5]. 

For high enough temperatures, when 27rT 3> u> s f, 
spin dynamics can be neglected [4] and electron 
interaction with spin (pseudogap) fluctuations re- 
duces to scattering by appropriately defined static 



Gaussian random field. For such model we can 
further simplify the form of effective interaction 
(correlator of the random field) [4,5], allowing the 
complete summation of Feynman perturbation se- 
ries, leading to the following recurrence procedure 
determining single - electron Green's function as 
G(e n p) = G fc=0 (e„p) from: 



G k (e n p) = 



1 



- £fc(p) + ikv k K - E fc (e„p) 



E fc ( £ „p) = W 2 s(k + l)G k+1 (e nP ) 



(1) 



(2) 



where k = £ 1 — is the inverse correlation length 
of the pseudogap fluctuations, e n = 2irT(n + 1/2), 
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where v(p) = - is the velocity of quasiparticle 
with the spectrum £ p , which is taken in the stan- 
dard form [4] : 

£ p = — 2t(cosp x a + cospj,a) — At cosp x acosp y a — /i 

(5) 

where t — is transfer integral between nearest, while 
t' - between second nearest neighbors on the square 
lattice, fi — is the chemical potential. Parameter W 
with dimension of energy determines the effective 
width of the pseudogap. It is clear that both W 
and correlation length £ are (within our semiphe- 
nomenological approach) some functions of carrier 
concentration (and probably temperature) to be de- 
termined from the experiment [4,5]. 

The value of s(k) is determined by combinatorics 
of Feynman diagrams and for the case of Heisenberg 
spin (SDW) fluctuations [4] is equal to: 
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fc+2 



for odd k 



3 for even k 



(6) 



The limits of applicability of our approximations 
were discussed in detail in Refs. [4,5]. 

Remarkable advantage of our model is the pos- 
sibility of complete summation of all Feynman dia- 
grams (including those with intersecting interaction 
lines) also for the vertex parts, determining the re- 
sponse functions to an arbitrary external pertur- 
bation. Detailed enough discussion of the vertex 
parts was presented in Ref. [8]. Here we just write 
down the appropriate recurrence relations appear- 
ing after the appropriate analysis for "triangular" 
vertices in Cooper channel, analogous to those de- 
rived in Ref. [6] and describing the response to an 
arbitrary fluctuation of superconducting order pa- 
rameter (energy gap) A q e(p), where the symme- 
try factor determining the type (symmetry) of pair- 
ing is taken here for the case of d-wave pairing as 
e(p) = cosp^a — cosp^a, and we assume the usual 
case of singlet pairing. The vertex of interest to us 
can be written as: 

r(e„, -e n , P, -P + q) = r p (e„, -e n , q)e(p) (7) 

Then r p (e n ,— e n ,q) is determined by the recur- 
rence procedure of the following form: 

r p fc_i(e n ,-e n ,q) = 1 - W 2 r(k)G k (e n , p + q) x 
xG k (-e n , p)S fe (e„, p, q)T pk (e n , -e n , q) (8) 

where 



S fe (e„,p,q) = { 1 



+ 



2iknv k 



G k 1 (s n , p + q) - G k 1 (-e n , p) - 2iknv k 



(9) 



"Physical" vertex is given by T pk= o(e n , ~e n , q). 
Additional combinatorial factor r(k) for the case of 
Heisenberg spin (SDW) fluctuations [4] is given by: 



r(fc) = 



k 

fc+2 
9 



for even k 
for odd k 



(10) 



The choice of the sign before the factor of W 2 in the 
right hand side of Eq. (8) depends on the symmetry 
of superconducting order parameter and the nature 
of pseudogap fluctuations [6] , here our choice again 
corresponds to the case of Heisenberg spin (SDW) 
fluctuations [10]. 

Scattering by normal (non magnetic) impurities 
can easily be accounted for by self - consistent Born 
approximation writing down "Dyson equation" for 



single - electron Green's function generalizing Eq. 
(1) and shown diagrammatically in Fig. 1 (a), 
where we have introduced: 



G fc(e n p) 



1 



(11) 



i£ n - 6c (p) + ikv k K 
so that instead of (2) we get: 

S fc (e„p) = P U 2 J2 G(enP) + W 2 s(k + l)Gfc+i(e n p) 
p 

(12) 

where p — is concentration of point - like impurities 
with potential U, and impurity contribution to self 
- energy contains the fully dressed Green's function 
G(e„p), to be determined self - consistently. 
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FIG. 1. Recurrence relations for the Green's function 
(a) and "triangular" vertex (b) with the account of scat- 
tering by random impurities. 

In comparison with impurity free case, we have 
an obvious substitution (renormalization) : 



P U 2 IrnG{e n p) 



(13) 



Dropping full self - consistency of self - energy over 
impurity scattering we obtain just the usual substi- 
tution: 
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+ 7osigns n (14) 

where 70 = irpU 2 No(Q) - is standard Born impu- 
rity scattering rate (iVo(O) - free electron density of 
states at the Fermi level). 

For a "triangular" vertex with the account of im- 
purity scattering we obtain a recurrence relation of 
general form shown graphically in Fig.l (b). How- 
ever for the case of the vertex, describing interaction 
with a fluctuation of superconducting order parame- 
ter with d-wave symmetry this equation is simplified 
considerably as the second term in the right hand 
side of Fig. 1 (b) is in fact zero due to ^ p e(p) = 
(cf. similar situation discussed in Ref. [9]). Then 
our recursion relation for the vertex part reduces to 
(8), where Gfe(±e n p) are determined by Eqs. (11), 
(12), i.e. are just Green's function "dressed" by im- 
purity scattering and defined diagrammatically in 
Fig. 1 (a). 

Superconducting transition temperature is deter- 
mined by the usual equation for Cooper instability 
of the normal phase: 

l-y X (0;T) = (15) 

where the generalized Cooper channel susceptibility 
is given by: 

x(q;T) = -T£5> 2 (p)x 

n p 

xG(e„p + q)G(-e„, -p)r p (e„, -e„, q) (16) 

Pairing interaction constant V is, as usual in BCS 
- like approach, assumed to be non zero in some 
layer of the width of 2u> c around the Fermi level 
and determines the "bare" transition temperature 
T c q in the absence of pseudogap fluctuations by the 
standard BCS equation [6]. Then, choosing rather 
arbitrarily e.g. uj c = OAt and T c0 = O.Oli we can 
easily find the appropriate value of V corresponding 
to this value of T c q. Actually this can be easily done 
for any other choice of parameters of our model. 

Knowledge of Green's functions and vertices al- 
lows us to derive Ginzburg - Landau expansion in 
a standard way [6] . Then, using microscopic values 
of Ginzburg - Landau coefficients we can determine 
all characteristics of a superconductor close to T c 
(e.g. coherence length, penetration depth, specific 
heat discontinuity etc.). Appropriate results will be 
given in a more detailed publication [10]. 

Numerical calculations were performed for differ- 
ent values of t, t', and ji of the "bare" electronic 
spectrum and also for different values of W and in- 
verse correlation length na. 

In Fig. 2 we show the dependence of super- 
conducting critical temperature T c on the effective 



width of the pseudogap W for different values of 
impurity scattering rate. It is seen that pseudo- 
gap fluctuations lead to significant suppression of 
superconductivity and that in case of finite im- 
purity scattering we always obtain some "critical" 
value of W, corresponding to complete disappear- 
ance of T c . This T c suppression is obviously due 
to partial "dielectrization" of electronic spectrum 
around the "hot spots" on the Fermi surface [4,5]. 
T c dependence on the value of correlation length 
of short - range order (pseudogap) fluctuations is 
much slower, in all cases the growth of £ (dimin- 
ishing n) enhances the effect of pseudogap fluctua- 
tions (suppression of T c ). We drop the appropriate 
results to spare space. 

Now we can move to an attempt to model a 
typical phase diagram of a cuprate superconductor 
within "hot spots" model of the pseudogap state. 
First attempt of such modelling in an oversimpli- 
fied ( "toy" ) version of our model was undertaken in 
Ref. [11] . The main idea is to identify our parameter 
W with experimentally observable effective width 
of the pseudogap (temperature of crossover to the 
pseudogap region of the phase diagram) E g « T* , 
which is determined from numerous experiments 
[1-3]. It is well known that the value of this pa- 
rameter actually drops, more or less linearly with 
concentration of doping impurity (carrier concen- 
tration), from values of the order of 10 3 K, becoming 
zero at some critical concentration x c « 0.19. .0.22, 
slightly greater than the "optimal" concentration 
x Q w 0.15. .0.17 [1,12]. Accordingly we just assume 
similar concentration dependence of our effective 
pseudogap width W(x). In this sense we can con- 
sider our W{x) as determined from the experiment. 
Then the only "fitting" parameter of the model is 
concentration dependence of the "bare" supercon- 
ducting transition temperature T c0 (x), which would 
have been existent in the absence of pseudogap fluc- 
tuations. Unfortunately, as already has been noted 
in Ref. [6], this dependence of T c0 (x) is, in gen- 
eral case, unknown and is not determined from any 
known experiment. 
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W/Tco 

FIG. 2. T c dependence on the effective width of the 
pseudogap W for the case of d-wave pairing and different 
values of impurity scattering rate 70/Tco: 0-1; 0.18 - 
2; 0.64 - 3. Inverse correlation length na=0.2. 

In the framework of our BCS - like approach any 
significant concentration dependence of the "bare" 
T c q seems to be unrealistic 1 . So we just assume 
there is no dependence of T c q on x at all, but take 
into account the fact, that introduction of any dop- 
ing impurity obviously leads to appearance of some 
impurity scattering (internal disorder), which can 
be described by the appropriate linear dependence 
of 7(2:). Assume that it is this growth of disorder 
is what leads to the complete suppression of d-wave 
pairing (according to the well known Abrikosov - 
Gorkov dependence, cf. e.g. Ref. [9]) say at x = 0.3. 
Then our approach allows to calculate concentra- 
tion dependence T c (x) for all values of x. Results 
of such calculation with parameters more or less 
appropriate for La 2 - x Sr x Cu04, in case of Heisen- 
berg (SDW) pseudogap fluctuations, and with the 
account of this very simple model of impurity scat- 
tering, are shown in Fig. 3. Actual values of pa- 
rameters of the model used in these calculations are 
shown at the same figure. "Experimental" data for 
T c (x) shown in this figure with "diamonds" are ob- 
tained from an empirical relation [12,13]: 



T c (x) 



1 - 82.6(2; - x y 



(17) 



served concentration dependence of T c for a number 
of cuprates. We can see that in the whole under- 
doped region our model gives practically ideal fit 
to "experimental" data, assuming quite reasonable 
values of W(x) 2 . 

It is interesting to analyze dependence of super- 
conducting transition temperature T c under addi- 
tional disordering of the system at different compo- 
sitions (concentration of carriers) . Such disordering 
was studied in a number of experiments, e.g. in- 
troducing disorder by neutron irradiation [14] or by 
chemical substitutions [15], in our model it can be 
simulated by additional impurity scattering, char- 
acterized by an extra parameter 70, which is just 
added to our parameter of an internal disorder 7(2;). 
Results of this type of calculations for two values of 
this additional scattering parameter are also shown 
in Fig. 3. It is seen that in complete accordance 
with experiments of Ref. [15], introduction of addi- 
tional "impurities" (disorder) leads to rather fast re- 
duction of superconducting region on the phase dia- 
gram. Also we note that in accordance with exper- 
iments [14,15], suppression of superconductivity by 
disorder in the underdoped region of the phase dia- 
gram (pseudogap state) is significantly faster, than 
at optimal dopings. 

This work was supported in part by the grant A'" 
02-02-16031 from the RFBR as well as by the RAS 
Programs "Quantum macrophysics" and "Strongly 
correlated electrons in semiconductors, metals, su- 
perconductors and magnetic materials" . 



2 Note that in case of Ising like pseudogap fluctuations 
we need rather unrealistic values of W(x = 0), order 
a magnitude larger, than in Heisenberg case. However, 
CDW like fluctuations can also give rather good fit [10] . 



Tc {x — x ) 

which gives rather good fit to experimentally ob- 



x This can be only due to some relatively smooth con- 
centration dependence of the density of states at the 
Fermi level. 



Journal of Physics and Chemistry of Solids 



250 



200 



150 



t*4 



100 



50 



Win(x=0)=580K; Tco=70K 
W(x) t'/t=-0,25, n/t^0,8, Ka=0,2; 
Tmax(x=0, 1 6)=39K 



1 - Y =0, 1 5Tco 

2 - y ( =0,25Tco 




o 

0,00 0,05 0,10 0,15 0,20 0,25 0,30 



FIG. 3. Model phase diagram for La,2- x Sr x Cu04 - 
type system in the case of concentration independent 
"bare" T c o with "internal" disorder j(x) linear in x. 
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